THE NUMBER OF NOWHERE-ZERO FLOWS 
ON GRAPHS AND SIGNED GRAPHS 



MATTHIAS BECK AND THOMAS ZASLAVSKY 



Abstract. A nowhere-zero fc-flow on a graph T is a mapping from the edges of T to the 
set {±1, ±2, . . . , ±(fc — 1)} C Z such that, in any fixed orientation of T, at each node the 
sum of the labels over the edges pointing towards the node equals the sum over the edges 
pointing away from the node. We show that the existence of an integral flow polynomial 
that counts nowhere-zero fc-flows on a graph, due to Kochol, is a consequence of a general 
theory of inside-out polytopes. The same holds for flows on signed graphs. We develop these 
theories, as well as the related counting theory of nowhere-zero flows on a signed graph with 
values in an abelian group of odd order. Our results are of two kinds: polynomiality or 
quasipolynomiality of the flow counting functions, and reciprocity laws that interpret the 
evaluations of the flow polynomials at negative integers in terms of the combinatorics of the 
graph. 
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1. Introduction 

A nowhere-zero flow on a graph F = (V,E), with values in an abelian group A, is a 
mapping x : E — > A such that, for every node v £ V, 

£ x(e) = x ^ ( L1 ) 

h(e)=v t(e)=v 

and which never takes the value 0. Here h(e) and t(e) are respectively the head and tail of 
the edge e in a (fixed) orientation of T. (In a certain sense, described below, x is independent 
of the chosen orientation.) A nowhere-zero k-flow is an integral flow (i.e., A — Z) whose 
absolute values are in [k — 1] := {1, 2, . . . , k — 1}. Nowhere-zero flows are nicely surveyed in 
PH] and P2]. 

It has long been known that the number of nowhere-zero flows with values in a finite 
abelian group of order k is a polynomial function of k (Tutte [22] )• Recently Kochol ^2] 
discovered that the number of nowhere-zero /c-flows is also a polynomial in k, although not 
the same polynomial. Here we show that this fact is a consequence of a general theory of 
counting lattice points in inside-out polytopes |3]. Furthermore, we extend Kochol's theorem 
in two ways: by a reciprocity law that combinatorially interprets negative arguments, and 
to signed graphs (in which each edge is positive or negative), where the polynomial becomes 
a quasipolynomial of period two: that is, a pair of polynomials, one for odd values of k 
and the other for even k ( Theorem I4.5|) : and we partially extend to signed graphs Tutte's 
counting theory for nowhere-zero flows in abelian groups (Theorem 14. Still further, the 
concept of reciprocity in lattice-point counting leads us to a geometrical interpretation of the 
number of totally cyclic orientations that are compatible with a given fc-flow ( Theorems 13.11 
and 14 .51 and Corollary I4.4|) . a fact which parallels Stanley's theorem [201 that the chromatic 
polynomial of T evaluated at negative integers counts acyclic orientations compatible with 
a node-labelling of T. 

2. The method of inside-out polytopes 

The theory of inside-out polytopes was motivated by the problem of counting the integral 
points of a rational convex polytope (the convex hull of finitely many rational points in M. d ) 
that do not lie in any of the members of a particular rational hyperplane arrangement. A 
(homogeneous, real) hyperplane arrangement is a finite set of homogeneous hyperplanes in 
M. d (that is, hyperplanes that contain the origin); it is rational if each of its hyperplanes is 
spanned by the rational points it contains or, alternatively, if each hyperplane has a rational 
normal vector. Suppose we are given a rational convex polytope P spanning M. d and a 
rational hyperplane arrangement IK. Then (P, Ji) is a rational inside-out polytope. More 
generally, P and "K may lie in a rational (and homogeneous) subspace Z that is spanned by 
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P. We introduced inside-out polytopes recently jl]; one motivation was to explain graph- 
theoretic enumeration functions geometrically. We consider inside-out theory essential to 
understanding our results on integral flows. 

A region (more precisely, an open region) of is a connected component of M. d \ |J !K; its 
closure is a closed region. The arrangement induced by Ji in a subspace S of R d is 

% s := {H r\ S : H e'K.H ^ S}. 
The intersection lattice of "K is 

£(!K) := {H§ : § C IK}, 

ordered by reverse inclusion j2H]; its elements are the flats of K. £ is a geometric lattice 
with 6 = f"| = R d and 1 = (~) K. (For matroids and geometric lattices we refer to ^B] or 
|21j.) The Mobius function of L is the function /i : £ x L — > Z defined recursively by 

r o if p g 5, 

/i(P, 5) := \ 1 if R = S, 

{-E R <u<sKR,U) i£R<S. 

(Sources for the Mobius function are, inter alia, ^H] and [2Tj .) 

2.1. Ehrhart theory. We outline the classical Ehrhart theory of lattice-point enumeration 
in polytopes. (See [7j, j5TJ Section 4.6, pp. 235-241], or [3].) We have a rational convex 
polytope P that spans a subspace Z. P° denotes the relative interior of P. The denominator 
of P is the least common denominator of all the coordinates of vertices of P. We denote 
by volP the volume of P, normalized with respect to Z PI Z d ; that is, we take the volume 
of a fundamental domain of the integer lattice in Z to be 1. To explain this last, we note 
that Z fl Z d is linearly equivalent to Z dimZ C M dimZ ; a fundamental domain is a domain in 
Z that corresponds to the unit hypercube [0,l] dimZ C M dimZ , under some invertible linear 
transformation that carries ZnZ d to Z dimZ . (See [2H pp. 238-239] or [3 Section 5.4] for 
more detail.) When Z = M. d this is the ordinary volume. 

A quasipolynomial is a function Q(t) = X^o Ci W^ defined on Z with coefficients c, that 
are periodic functions of t. Then Q is a polynomial on each residue class i modulo 
some integer, called the period of Q; these polynomials are the constituents of Q. (For an 
introduction to quasipolynomials see, e.g., j2U Chapter 4] or 

The subject of Ehrhart theory is the (closed) Ehrhart counting function 

E P (t) := #(tP D Z d ) 

and the open Ehrhart counting function, Ep°(t). Ehrhart's theorem is that Ep and i?p° are 
quasipolynomials with leading term (volP)t dimP and with periods that divide the denomi- 
nator of P. It follows that one can define both counting functions for negative integers; the 
Ehrhart-Macdonald Reciprocity Theorem is that Epo(t) = (—l) dimP Ep(—t). 

2.2. Inside-out Ehrhart theory. An open region of (P, "K) is a nonempty intersection 
with P° of an open region of Ji (thus it is full-dimensional in the span of P). A closed 
region of (P, !K) is the closure of an open region. A vertex of (P, IK) is a vertex of any such 
region. The denominator of (P, !K) is the least common denominator of all coordinates of 
all vertices. 
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The fundamental counting functions associated with (P, IK) are two quasipolynomials: the 
(closed) Ehrhart quasipolynomial, 



Ep,x{t) '■= ^2 m pMx), 
where the multiplicity mp^{x) of x G M. d with respect to IK and P is defined through 



the number of closed regions of (P, IK) that contain x, if x G P, 
0, if x <£ P, 



and the open Ehrhart quasipolynomial, 

E [P ,Ky(t) : = # (t-W n [p° \ ( |J IK) 

(This is simplified from the notation E po ^ of 4J. The reader may think of (P, IK) as the 
set P° \ (IJIK), called the relative interior of (P, IK), but that is not necessary in order to 
read this paper.) If P spans a subspace Z, then IK also lies in Z and the multiplicity of x is 
defined to be if x ^ Z. The names of our counting functions are justified by the fact that 
in the absence of IK we recover classical Ehrhart theory, and by one of the main results in 

® 

Theorem 2.1 ( 4, Theorem 4.1]). If (P, IK) is a closed, full- dimensional, rational inside-out 
polytope in Z C M. d , then Ep^(t) and E^p^o(t) are quasipolynomials in t, with period equal 
to a divisor of the denominator of(P,'K), with leading term (volP)t dimP , and with constant 
term E P ^(0) equal to the number of regions of (P, IK). Furthermore, 

E {P ^ r (t) = (-l) dimP Ep^(-t). (2.1) 

In particular, if (P, IK) is integral then Ep^ and Em^y are polynomials. The proof, 
though more general (and arrived at independently), is similar to Kochol's proof of Theorem 
13.1( a): we sum the Ehrhart quasipolynomials of the pieces into which IK dissects P. 

For the second theorem of |4| we will use here, we need the notion of transversality: IK 
is called transverse to P if every flat U G £(IK) that intersects P also intersects P°, and P 
does not lie in any of the hyperplanes of IK. 

Theorem 2.2 ( 4, Theorem 4.2]). If P and IK are as in Theorem \2.1\. then 

E {P M Y {t)= v{6,U)E P ° nU {t), (2.2) 



and if IK is transverse to P, 



E Pt x(t)= IM0,P)|P pnu(t). (2.3) 



Transversality is always satisfied in the applications here. 
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2.3. Matrix matroids. We shall want a general lemma about matroids of hyperplane ar- 
rangements induced by coordinate arrangements. We start with the hyperplane arrangement 
IK m , consisting of the coordinate hyperplanes in F m for some field F, and we take any sub- 
space S. Then S induces an arrangement 'K^ in S. (We are interested in the reals, but the 
lemma is valid for any field.) 

Any homogeneous hyperplane arrangement IK has a matroid M(!K), whose ground set 
is the set of hyperplanes and whose rank function is rkS = codim(f'lS) for § C "K. This 
matroid is simply the linear dependence matroid of the normal vectors of the hyperplanes. 
The column matroid of a matrix A, M(A), is the matroid of linear dependence of its columns; 
to keep the notation correct we take the ground set to be the set of indices of columns. The 
chain-group matroid of a subspace S C F m is the matroid N(S) on [m] whose circuits are 
the minimal nonempty supports of vectors in S. Lat M denotes the lattice of closed sets of 
a matroid M. Thus LatM(IK) ^ £(J£). 

We refer in the following lemma to orientations of oriented matroids. The reader who 
is not familiar with oriented matroids should think of orientations of a graph or (later, in 
Section Hj) of a signed graph. 

Lemma 2.3. Let A be an n x m matrix with entries in a field F , let Ji m = {H e : e G [m]} 
be the arrangement of coordinate hyperplanes in F m , and let U = Row A, the row space, and 
Z = Nul A, the null space. 

(a) The mapping e i— > H e D U is a matroid isomorphism from M(A) to M(IK^). Also, 
e i— > H e H Z is a matroid isomorphism of iV(Row A) with M^K^). 

(b) The mapping 

S i-> E V (S) := {e G [m] :H e DS} 

is the isomorphism of £(J{^) with Lat M(A) induced by the first mapping in (a). 
The mapping 

S i-> E Z (S) := {e G [m] : H e D S} 

is the isomorphism o/£(J{f) with Lat iV(Rowv4) induced by the second mapping in 
(a). 

(c) If F is an ordered field, then the regions of correspond bijectively to the acyclic 
orientations of the oriented matroid of columns of A, and those ofK^ correspond to 
the totally cyclic orientations. 

Proof, (a) Each coordinate hyperplane H e equals bjr for a basis vector b e along the e-axis. 
Therefore, H e D U is the orthogonal complement in U of a e , the orthogonal projection of 
b e into U. Let B := {b e : e G [m]}. We may take A to be the matrix whose columns are 
the vectors a e , since it has the same row space as the original matrix A. By linear duality, 
M(A) = M(J{^) under the correspondence e t— > H e fl U. 

We treat "K^ by taking A*, a matrix whose row space is Nul A, and applying the first part 
to Row A*. N(RowA) is dual to M(A); hence it is M(A*) because Row A and Nul A are 
dual chain-groups. (See, e.g., Tutte [22, Chapter VIII] on primitive chain-groups.) 

(b) Obvious from (a) and the definitions. 

(c) Write 'M.(A) for the oriented matroid. The proof is as in part (a), but relying on 
the dual relationship between M(A) and M(A*) [31 3 Section 6.3(c)]. Our M(A) is Ziegler's 
M(Nul A), so our M(A*) = M(Row A), which is dual to M(NulA) because the row and null 
spaces are orthogonal complements. ■ 
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3. Flows on graphs 



A flow on a graph F with values in an abelian group A, called an A-flow, is a function 
x : E — > A which satisfies for every node v G V (so it is like a nowhere- zero flow, but the 
flow value zero is allowed). This definition requires that the edges be oriented in a fixed way. 
The orientation is arbitrary; it is an artifact of notation, and to overcome this artificiality 
we define, for an oriented edge e, e _1 to be the same edge in the opposite orientation and 
^(e" 1 ) := —x(e). With this law for flows, the validity of Equation II. II is independent of the 
choice of the orientation of T. 

Tutte [22 Section 6] proved that the number of nowhere- zero A-flows on F is a polynomial 
in \ A\, independent of the actual group. We shall write (pr for this polynomial and call it the 
(strict) modular flow polynomial of F. (Usually (p? is called just the "flow polynomial" but 
we need to distinguish it from other flow polynomials.) As Tutte showed, the modular flow 
polynomial <pr(k) is the evaluation (— l)^ r )tr(0, 1 — fc) of the Tutte polynomial of F. (£(T) is 
the cyclomatic number \E\ — \V\ +c(F), where c(F) is the number of connected components.) 
We should mention that the total number of A- flows is the simple polynomial |A|^ r \ Tutte 
further proved (in [2U pp. 83-84], based on |2*3*1 Theorem VI]) that a nowhere-zero Z^-flow 
exists if and only if there is a nowhere-zero fc-flow, a k-flow being an integer-valued flow 
whose values all satisfy \x(e)\ < k. However, the actual number of nowhere-zero fc-flows for 
k > 0, which we write as <pr(k), does not equal the number of nowhere-zero Z fc -flows and 
indeed was never known to be a polynomial until the recent work of Kochol ^2J- Kochol 
employed standard Ehrhart theory combined with a special construction to prove this. We 
shall show that Kochol's theorem is a natural consequence of inside-out polytope theory and 
can be extended to a reciprocity theorem that interprets <fr{k) at negative arguments. 

A circle in a graph is a 2-regular connected subgraph. A cycle (or directed cycle) is a 
circle in which the edges are oriented in a consistent direction. An orientation of F is acyclic 
if it has no cycles and totally cyclic if every edge lies in a cycle. We call a totally cyclic 
orientation r and a flow x compatible if x > when it is expressed in terms of r. Taking 
the standpoint of the flow x, the nonzero edge set suppx has a preferred orientation, the 
one in which x > (we call this t(x)\ note that it orients only suppx) and the zero edges 
are free to take up any orientation that makes F totally cyclic. An isthmus of a graph is 
an edge whose deletion increases the number of connected components. There is no totally 
cyclic orientation if F has an isthmus. 

The real cycle space Z is defined in R E by Equation (jl.lj) . To this space Z we associate 
the polytope and arrangement 

P:=Zn[-l,l] E , X;=(X E ) Z , 

where "He is the arrangement of coordinate hyperplanes in M. E . A (k + l)-flow is then 
precisely a point x G Z n Z d such that tX G P and a nowhere-zero A;-flow is just a point 
x G Z n Z d such that |igP \ (|J IK). Consequently, 

<p° r (k + 1) = E P (k) (3.1) 

and 

(fr(k) — E(p 7 H) (k). (3.2) 

Theorem 3.1. Given: a graph F and its real cycle space Z . 

(a) (Kochol [T2~] ) <fir{k) is a polynomial function of k for k = 1, 2, 3, . . . . It has leading 

term (vol P)k^ r ^ ifF has no isthmi; otherwise it is identically zero. 
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(b) Furthermore, (—l)^ T ^<fr(~k) for k > equals the number of (k + 1) -flows counted 
with multiplicity equal to the number of totally cyclic orientations of F that are com- 
patible with the flow. 

(c) In particular, the constant term <^r(0) equals the number of totally cyclic orientations 
ofY, which equals (— l)^ r )<^ r (— 1). 

(d) Finally, the total number of k-flows for k > 0, nowhere-zero or not, is a polynomial 
Pr(k) satisfying ^(k) = (— l)^ r Vr(l — k), whose leading term is the same as that 
of (fr(k) and whose constant term is (—1)^ T \ 

We call ifr the (strict) integral flow polynomial of T and <p®{k) the weak integral flow 
polynomial. <^r(^) already gives rise to a number of interesting computational problems, as 
discussed, for example, in [2]. 

The reciprocity formula (p®(k) = (— l)^ r )<yjp(l — k) leads to the conclusion that tp® has no 
integral zeros. For any k > there is an integral flow, the all- zero flow. Reciprocity takes 
care of k < 0. 

Proof. For (a) we apply Theorem 12.11 in Z. We call upon the total unimodularity of the 
matrix of the cycle equations (jl.lj) to deduce that P is a convex hull of integer lattice points. 
Thus, Emjq° is a polynomial, and as we saw in (|3.2|) . it equals ipr- 

Since Ep t ^(k) counts pairs (x, R) where x G Z d H P and R is a closed region of 0~C that 
contains x, part (b) follows if we show that the regions of "K correspond to the totally cyclic 
orientations of V and a region of "K whose closure contains a chosen point x G Z H Z d 
corresponds to a totally cyclic orientation that is compatible with x. The first statement 
was demonstrated by Greene and Zaslavsky in jHJ Section 8], based on the obvious bijection 
(given the fixed orientation of T) between orthants of M. E and orientations of T. The second 
is then obvious. 

Thus the constant term is the number of totally cyclic orientations. The fact that this 
equals tr(0,2) is a theorem originally due to Las Vergnas (see [131 Proposition 8.1] and [HJ 
remark after Theorem 1', p. 296]) and independently proved by Greene and Zaslavsky (Hj 
Corollary 8.2]. 

Part (d) is standard Ehrhart theory, because a fc-flow is simply a point x G Z d such that 
-r£ G P°. That is, <^r(^) = Epo(k) = (—lY^Ep^—k) by Ehrhart reciprocity. The constant 
term of Ep(—k) is 1, the Euler characteristic of P. It is easy to see that Epo(k) = Ep(k — 1) 
for k > 0. Consequently, 

<p${k) = (-l) m E P (-k) = (-1)^E P °(1 -k) = (-lFVr(l - k) 

if k is a positive integer, whence for all k. ■ 

Problem 3.2. Find a formula for, or a combinatorial interpretation of, the leading coefficient 
vol P of the integral flow polynomials. 

Problem 3.3. Is there a combinatorial interpretation of <^r( — k) for k > 2? 

Example 3.4 (Small graphs). We calculated the integral flow polynomials of some small 
graphs by counting integral /c-flows on a computer for 1 < k < £ (r) + 2 and interpolating to 
get the polynomial. The graphs were mi^, the graph of m parallel links, for m = 3,4, 5, 6, 
and K^. We state our results along with the modular flow polynomials for comparison; the 
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latter are (fr(k) = k^ and <pr(k) = Xr*(k)/k, T* being the planar dual graph. First, 3K 2 : 

^°{k) = k 2 , (p(k) — {k — l)(k — 2), 

ip\k) = 3k 2 -3k + 1, <p{k) = 3(k - l)(Jfe - 2). 

Next, AK 2 : 

f(k) = k 3 , Cp{k) = (k — l){k 2 -3k + 3), 

0/r . (2k- l)(8k 2 -8k + 3) ,,, 2(k - l)(8k 2 - 22k + 21) 

V C«J = o ' <PW = 5 • 



Next, 5K 2 : 



(p°(k) = k\ (p(k) = (k - 1)(A; 3 - Ak 2 + Qk - 4), 
115A; 4 - 230A; 3 + 185A; 2 - 70A; + 12 



12 

5(Jfe - l)(Jfe - 2)(23A; 2 - 41Jfe + 36) 
12 ' 



Next, 6K 2 : 



(p°(k) = k 5 , (p(k) — (k — l)(k 4 - 5k 3 + 10k 2 - 10k + 5), 



<p(k) 



2{2k - 1)(44A; 4 - 88A; 3 + Ilk 2 - 27k + 5) 

io ' 

(k - l)(176/fc 4 - 839A; 3 + 1571A; 2 - 1404A; + 620) 
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Finally, K±: 

^°{k) = k 3 , Cp{k) — (k — l)(k - 2)(k - 3), 

(p°(k) = (2k - l)(2k 2 -2k + 1), (p(k) = 4(fc - l){k - 2){k - 3). 

Problem 3.5. Is there any general reason why in some of these examples {3K 2 and K4) both 
of the integral flow polynomials have integral coefficients and the integral nowhere-zero flow 
polynomial is a multiple of the modular polynomial? 

The totally cyclic orientations that are compatible with a flow x are the totally cyclic 
extensions of r{x) to V. We get such an extension by orienting the contraction V/ suppx in 
a totally cyclic manner. The number of such orientations is given by Theorem I3.1f c). Thus 
we have the following version of Theorem I3.1f b): 

Corollary 3.6. |<^r(-fc)l = J2 xe pnk-^ |£r/suw>*(-l)|- ■ 

The multiplicity of an integral point x with respect to a rational hyperplane arrangement 
IK is the number of closed regions of !K that contain x. As we pointed out in proving Theorem 
13.1( b). the number of closed regions of !Kf that contain a flow x is equal to the number of 
totally cyclic extensions to T of r{x). Hence, the multiplicity of the integral flow x with 
respect to CKf is another interpretation of the quantity (— l)^ ( ' r ^ suppa '' ) <^r/suppx(— !)• 

Other formulas arise from the intersection expansions of Theorem 12.21 but as we need its 
Mobius function, first we have to find the lattice £(!H) explicitly. We do so in the more 
general context of signed graphs. 
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4. Flows on signed graphs 



The best way to understand the cycle equations (jl.lj) is in terms of the incidence matrix, 
which we expound in the general context of signed or bidirected graphs. 

Formally, a signed graph £ = (T, a) consists of a graph T and a function a from the set of 
links and loops of T to {+, — }. (A link has two distinct endpoints; a loop has two coinciding 
endpoints. In signed and bidirected graph theory it is convenient to have two more kinds of 
edges: a half edge has one endpoint and a loose edge has no endpoints; neither of these has a 
sign.) If T C E, then T\T or E|T denotes the spanning subgraph whose edge set is T. Each 
circle has a sign, which is the product of the signs of its edges. A subgraph or edge set is 
called balanced if it contains no halfedges and every circle in it has positive sign. (See jU] for 
the origin of signed graphs and balance; for the general theory of signed graphs see [27|.) 

The bias matroid (or signed- graphic matroid) of S j2Z| Section 5 and Erratum], written 
G(S), can be defined by its rank function, 

r(T) = \V\- b(E\T) for an edge set T, 

where 6(S|T) is the number of components of the subgraph E|T that are balanced subgraphs, 
ignoring any loose edges. The circuits of G(E) are of three kinds: a positive circle, a pair of 
negative circles that have a single common node, or a pair of node-disjoint negative circles 
together with a minimal connecting path; here one or both negative circles may be replaced 
by halfedges. A coloop of G(Y0) is an edge e whose deletion makes an unbalanced component 
balanced; or which is an isthmus connecting two components of S \ e of which at least one 
is balanced (all this by |27J Theorem 5.1 as corrected]). We define the cyclomatic number of 
£ to be \E\ — \V\ + b(E). This is the rank of the dual G _L (E) of the bias matroid. 

A graph is bidirected when each end of each edge is independently oriented. We express 
the bidirection by means of an incidence function rj defined on the edge ends: the function 
is +1 if the arrow on that end points into the incident node, and —1 otherwise. For an edge 
end e, let e(e) denote the edge and v(e) the node incident to e. We define 

r](v, e) := {v( £ ) '■ v ( £ ) = v an d e(e) = e}. 

Thus, r](v, e) = if v and e are not incident or e is a loose edge or positive loop. The 

bidirected incidence matrix is 

E(T, V ) := {r)(v,e)) VxE . 

A bidirection of a graph is really an orientation of a signed graph. A link or loop e with 
ends e\ and e 2 has sign 

(7(e) : = - v ( £l )ri(e 2 ). 

In plain language, if the two arrows on e point in the same direction, then e is positive, 
but if they conflict, e is negative. We call 77 an orientation of the signed graph S. (See 
[2H] . This notion corresponds to the ordinary notion of graph orientation if we identify an 
unsigned graph F with the all-positive graph +F) We call an (oriented) incidence matrix 
of £ any incidence matrix of an orientation of S, writing it H(E) (and H(r) := H(+r)). 
The ambiguity arising from the unspecified orientation can usually be ignored, since the only 
effect of changing the orientation is to negate some columns. Recall that reorienting e to 
e _1 also has the effect of replacing x(e) by x(e~ l ) = —x(e) for every x G A E (A an abelian 
group). With these conventions we define a flow on E with values in A as any x G A E for 
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which 



H(E)x = 



(4.1) 



in other words, for which x G Nul H(£). This definition generalizes that of a flow on a graph; 
naturally, then, we generalize theorems about flows. A k-flow is an integral flow x for which 
every \x(e)\ < k, just as before. Nowhere-zero fc-flows on signed graphs were studied by 
Bouchet [5 and Khelladi [TT], searching for the smallest possible k. 

A cycle in an oriented signed graph is a circuit such that no node has all arrows pointing 
into the node or all arrows pointing out of the node. An orientation is acyclic if it has no 
cycles and totally cyclic if every edge belongs to a cycle. 

Suppose we have a bidirected graph. Switching a node v means changing rj to rj v defined 



The associated switched signed graph is denoted T, v . It is obtained from E by negating all 
links incident with v. 

4.1. Group- valued flows. We begin our treatment of signed graphs with the analog of the 
modular flow polynomial, since as far as we know it has not been published. The Tutte 
polynomial t^(x,y) of E is defined to be that of G(E). (It is not equal to that of the 
underlying graph unless E is balanced. See |Uj for the Tutte polynomial of a matroid.) 



Theorem 4.1. For each signed graph £ there is a polynomial <p^(k) such that the number 
of nowhere-zero flows on £ with values in a finite abelian group A of odd order is ip^(\A\). 
In fact, 



Proof. For our proof by induction on \E\ we need deletion and contraction of edges in signed 
graphs, as in [27] . Deletion is simply removing the edge with no other change. The contrac- 
tion E/e by an edge e is done differently for the four kinds of edge. A loose edge or positive 
loop is simply deleted. A halfedge or negative loop is deleted and its supporting node is also 
deleted, but any other edges incident to that node are retained; such an edge that was a link 
becomes a halfedge, but if it was a loop or halfedge it becomes a loose edge. A positive link 
is deleted and its endpoints are identified to a single node. To contract a negative link we 
switch an endpoint so the link becomes positive; then the link is deleted and the endpoints 
are identified. 

In the proof we fix A and an orientation 77 of £ and write (p-^{A) for the number of 
nowhere-zero A-flows on £. 

It is clear that if £ has components Si, ... , £ c , then <px(A) = CpsxiA) ■ ■ -<fs c (A). Analo- 
gously, ts = tsi • • 4 ts c . Therefore, we may assume £ is connected. 

Suppose e is a positive link of E, oriented from v to w, such that deleting e from E does 
not increase the number of balanced components; that is, e is not a coloop. Let 
nowhere-zero A-flow on E/e. We can define a flow on E by 



by 




x(f) = x"(f) iff^e, 
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Then H(£)sc = 0. If x(e) 7^ 0, x is a nowhere-zero flow on E; if x(e) = 0, a; is a nowhere- zero 
flow on E \ e. Conversely, any nowhere- zero flow on E \ e or E gives rise to one on E/e. 
Since the assumption about e guarantees that E, E/e, and E \ e all have the same number 
of balanced components, 

(-l)^MA) = (-lp /e W04) + (-l) e(EV W(A). (4.2) 

If e is a negative link, subject to the same hypothesis on balance, with endpoints v and 
i/, we switch v to make e positive. Because a flow on E remains a flow after switching, 
<Py?>{A) = (ps(A). Thus, (14. 2 j) is valid for a negative link. 

The analogous Tutte-polynomial formula, is = ts/e + ^s\e> is valid as long as e is not a 
coloop in G(E), but that is our balance assumption. In this case, therefore, the theorem is 
valid for E by induction. 

If on the other hand an edge e is a coloop, then ts(0, y) = and, by the proofs of Bouchet's 
Lemmas 2.4 and 2.5 also <^e(^4) — 0. (The oddness of A, so that 2a = implies a = in 
A, is needed at this point.) 

The preceding arguments reduce the theorem to the case of a signed graph with a single 
node v. The flow in a loose edge or positive loop can be any nonzero value in A, independent 
of all other flow values. Therefore, <fi:(A) = (\A\ — l) l (p-£, {A) where I is the number of 
positive loops and loose edges and E is E with all such edges deleted. Now let us assume 
E has one node and its edges are halfedges ei, . . . , e^ and negative loops fi, . . . , fj. We may 
assume they are oriented into v, so the total inflow at v is 

x(ei) + • • • + x( ei ) + 2x{h) + ••• + 2x{f j ) = 

by the cycle condition (j4.1j) . If % +j = 0, there is one such flow, in agreement with £g = 1. If 
i+j = 1, there is none (again we require \A\ odd) so <p% — 0, in agreement with t%(Q,y) = 0. 
If i + j > 1, we get a nowhere-zero A-flow by taking arbitrary nonzero values for x except 
on one edge, say x(e), which is determined by the rest through the cycle equations. The 
number of such choices in which x(e) 7^ is (px,(A); the number in which x(e) = is (ps\ e (A). 
Therefore, 

<fr{A) = - \A\y+i- 1 - <p* v {A) 

= (-lF /e) t S/e (0, 1 - \A\) - (-l)^ e) t sv (0, 
by the fact that G(E/e) has rank and induction on i + j = \E\, 

= (-1)^(0,1-14) 

by the fact that e is not a loop or coloop in G(E). Our theorem now follows by induction. ■ 

Problem 4.2. Is there any significance to or evaluated at even natural numbers? 

Theorem 14. II means there is a polynomial <p^(k), which we call the (strict) modular flow 
polynomial, such that for any odd positive number k, <fz{k) is the number of nowhere- 
zero flows on E with values in any fixed abelian group of order k. This is reminiscent 
of signed-graph coloring, where only odd values of A can be interpreted in the chromatic 
polynomial Xs(-^) |28| . For coloring, though, there is another polynomial which counts 
restricted colorings when evaluated at even integers; and the two polynomials are related (as 
we showed in |3]). We wonder whether there could be something similar with the modular 
flow polynomial, and whether flows and colorings might be connected through duality of 

signed graphs, analogously to the duality of colorings and flows on planar graphs. 
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Example 4.3. The number of A- flows depends on the group A when it has even order. 
Consider the signed graph consisting of two negative loops at a node v. Orient both loops 
into v. The cycle equation is 2x(e) + 2x(f) = 0. We compare two groups of order 4. 
If A is the Klein four-group Z2 x Z2, then every x : E — > A satisfies the cycle equations and 
we have 

<pl(A) = 16, MA) =9. 
HA — Z 4 , then the A-flows are (x(e), x(f)) e {0, 2} 2 U {1, 3} 2 and 

vl(A) = 8, MA) =5. 

For comparison, at odd values of k we have 

Mk) = k, Mk) = k-1. 

Corollary 4.4. The number of totally cyclic orientations o/E equals (—l)^ 2 ) </?£(— 1). 

Proof. The number of totally cyclic reorientations of an orientation of a matroid M is £m(0, 2) 
[T3*| IT4*] . Since cycles in an orientation of £ are the same as cycles in the corresponding 
orientation of G(£) [2E], the number of totally cyclic orientations of £ equals £g(£)(0,2) = 
(_1)€(=)^,(_1). ■ 

4.2. Integral fc-flows on signed graphs. Now it is time for integral flows. For k > let 

<pT,{k) '■= the number of nowhere-zero fc-flows on E. 

As with abelian-group flows, ^2 = if there is a coloop in G(E). Let 

<£>jj(A;) := the number of all fc-flows on E, 

also for k > 0. We take Z to be the real cycle space NulH(E) — this is the solution space of 
Equation (jl.lj) — and, just as with unsigned graphs, 

P := Zn[-l,lf and "K := ttf 

where 'He is the coordinate-hyperplane arrangement. As with ordinary graphs, a flow x and 
an orientation 77 are compatible if 2 > when expressed in terms of 77. 

Theorem 4.5. (oj For any signed graph E, ips(k) is a quasipolynomial function of k for 
k = 1, 2, 3, ... . Its period is 1 or 2, and is 1 if '£ is balanced. <pj;(k) has leading term 
(vol P)k^> if G(Y>) has no coloops; otherwise v?e(^) is identically zero. 

(b) Furthermore, (—lY^(fs(—k) for k > equals the number of {k + l)-flows counted 
with multiplicity equal to the number of compatible totally cyclic orientations o/E. 

(c) In particular, the constant term (/?s(0) equals the number of totally cyclic orientations 
o/E, which equals (— l)^( s ^ s (— 1). 

(d) Finally, <p^(k) is a quasipolynomial of period 1 or 2 (period lifEis balanced) whose 
leading term is the same as that o/^s(/c) and whose constant term is (— l)^ r ). Fur- 
thermore, M k ) = (-l) $(r) <Ps(l - k )- 

From (d) we see that, even when ip° has period 2, in a way it consists only of one polyno- 
mial. The odd constituent, v^odd; is determined by the even constituent through 

^dd(2j + i) = (-iFVLn(-2j)- 
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Problem 4.6. Is there a nontrivial unbalanced signed graph for which has period 1? (For 
the dual function, the chromatic quasipolynomial, there is not. See the second remark after 
[3J Corollary 5.9].) The graph of one node in Example 14.81 shows that trivial examples exist. 
As for there exists at least one such graph with more than one node, but still very small; 
see Example 14.81 

The mere existence of such examples is significant: it suggests that the duality between 
flows and colorings is imperfect and limited. 

Lemma 4.7. The vertices of (P, "K) are half integral. 

Proof. We may as well assume £ has a halfedge at every node; thus H(E) contains an 
identity matrix I n . A vertex is a solution of H(S)x = with | _E7| — rz coordinates of x set 
equal to fixed values in {0, 1, —1}. Let B be the edge set whose coordinates in x are left 
undetermined, let B c := E\ B, and write x = (xb,xb c ) T ■ Then x is the unique solution of 
H(S| J B)x B = -H(E|5 c )x S c. 

The remainder of the proof is based on work of Jon Lee. Lee proved that the null space 
NulH(S) is 2-regular Proposition 9.1] and that if A is a nonsingular square matrix for 
which Nul[J | A] is 2-regular, then A~ l b is half integral for every integral vector b (a special 
case of Proposition 6.1]; the definition of 2-regularity need not concern us). These facts 
applied to A = H(S|5) imply that the solution of H(E\B)x = b is half integral for any 
b e Z B . Apply this to b = -E(Z\B c )x B c ■ 

Proof of Theorem |^.<5| The proof is similar to that of Theorem 13.11 In (a) and (d), instead 
of total unimodularity we have Lemma f4. 71 to tell us that the denominator of (P, !K), hence 
the period of the Ehrhart quasipolynomials, divides 2. 

For (b) we need to show that the regions of 'K correspond to the totally cyclic orientations 
of S. The latter are the totally cyclic reorientations of the natural orientation of G(S), which 
is the oriented matroid of columns of H(S) Theorem 3.3]. Now we apply Lemma f2.3f b). 

For (c) we use Corollary 14.41 ■ 

Example 4.8. We calculated the flow polynomials and modular flow polynomials of some 
small signed graphs. First we treat the signed graph with two negative loops at one node, 
for comparison with the modular flow polynomials in Example 14.31 

<fl(k) = 2k - 1, y> E (Jfe) = 2(k - 1). 

Another signed graph, only slightly larger, has two nodes joined by positive and negative 
edges and at each node either a halfedge or a negative loop. We write iK^ for this graph 
if i nodes have halfedges and j other nodes have negative loops and <f(i,j) for the polynomials; 
the examples we calculated are where i + j = 2. The modular nowhere- zero flow polynomial 
is the characteristic polynomial of the dual matroid of Gi^kK^ ), which is the four-point 
line. The modular polynomials are 
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The integral polynomials are 

'2(k- l)(Ac — 3) if k is odd, 



^(0,2) (A;) = V(2,o)(k) 



2(k — 2) 2 if k is even 



and 



±(3fc 2 -2A; + l) if jfe is odd, 
±(3£; 2 -4A; + 2) if A; is even, 

|(3fc 2 - 12*; + 9) if k is odd, 
±(3£; 2 - Uk + 16) if A; is even. 
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Problem 4.9. We do not understand why ±K^' 2 ^ an d ±K^' ^ have the same flow quasipoly- 
nomials but ii^ 1 '^ does not, nor why (p® Q ^ has period one, nor why the odd constituent of 
<P(p,2) is an integral multiple of v?(m> 

4.3. Half integrality and the incidence matrix. In our geometric treatment of graph col- 
oring in |U Section 5] we noticed an important half-integrality property of the signed-graphic 
hyperplane arrangement CK[E], which consists of the dual hyperp lanes to all the columns in 
the incidence matrix of S. It is the next lemma, which involves a polytope P = [0, l] d and 
the affine arrangement J{"[£] which is !K[S] translated by the vector |(1, 1, . . . , 1) T . 

Lemma 4.10 (jU Lemma 5.7]). If S is a signed graph, (P, !K"[£]) has half-integral vertices. 

Despite their appearance, not only are Lemmas 14.71 and 14.101 similar in statement, both 
also involve the incidence matrix. A basic solution of a linear program H(S) T x < b, x > 
is the solution of H(S|5) T x = b for some set B of n linearly independent columns of H(£). 

Proposition 4.11. Every basic solution of H(S) T x < b, x > is half integral for any 
b E Z n . 

Proof that Lemma \4-lU\ o,nd Proposition \4-H\ are equivalent. We may assume that £ con- 
tains a negative loop at every node. Write a e for the column of the incidence matrix that 
belongs to e. The equations of the hyperplanes in U-CfE] £1X6 Q> e X — 

(1, 1, . . . , 1) ; this is 

or 1 if e is a loop or link but it is | if e is a halfedge at node Vi and in that case the equation 
simplifies to Xi — |. In the situation of Lemma 14.101 t ne vertices of (P, IK"[£]) are obtained 
by choosing from among the equations x^ — 0, x^ — 1 (from P) and Xi = |, ajx = or 1 
(from CK"[S]) a subset of n equations whose coefficient matrix is nonsingular. An equation 
Xi = or 1 has the form ajx = or 1 with e the negative loop at Vi, so we need not consider 
it separately. Thus, half integrality of vertices of ([0, l] n , [K"[E]) implies half integrality of 
the solution of H(S|i?) T x = b for any b e Z n . The converse is obvious. ■ 

Our attention was drawn to the relationship between Lemma 14.101 and half integrality 
(and Lee's work) by a recent manuscript of Appa and Kotnyek [l]. 



4.4. Nowhere-zero flows reduce to flows with Mobius complications. The last main 
result expresses the nowhere-zero integral flow polynomial in terms of the weak integral flow 
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polynomials of subgraphs. We begin with structural lemmas. As before, Z is the cycle space, 
P := [-1, 1} E n Z, and "K := ttf . For a flat S G tL{K) we define 

:= {e G £ : 5 C H e } = {e G E : ar(e) = for all x G 5}. 

This is the E Z (S) of Lemma l2~37 b). We see that E(S) C is the union of the suppports of the 
vectors in S. 

Lemma 4.12. The lattice of flats of"K is isomorphic to the lattice of closed sets of the dual 
of the bias matroid, G- L (S). The isomorphism is given by S i— > E(S). The corresponding 
matroid isomorphism G _L (S) = M{ { K) is given by e i— > H e D Z . 

Proof. This is an application of Lemma f2. HI The matrix is H(S). We know M(H(S)) is the 
bias matroid G(H) by fZT\ Theorem 8A.1], so G _L (E) is the chain-group matroid of Row H(S). 
The lemma applies because the real cycle space Z = NulH(S). ■ 

(Lemma 14. 121 holds good for the canonical hyperplane representation of any F*-gain graph 
<3>, for any field F [101 Section 2}. Denoting an incidence matrix by H($) and taking Z = 
NulH($), we have Af(J{|) = But we digress.) 

Lemma 4.13. (fs is identically zero if and only if Gill) has a coloop. 

Proof. This is an application of Bouchet's theorem on integral chain-group matroids jSJ 
Proposition 3.1]: the chain-group has a nowhere- zero chain if and only if the matroid has no 
coloop. In our case the chain-group is the group of integral flows, NulH(S). Its chain-group 
matroid is dual to that of RowH(S), which is dual to the column matroid of H(E), which is 
G(Yj). ■ 

Lemma 4.14. A flat S ofK can be represented as [Nul H(E|E(5) C )] x {0} E ( S \ 

Proof. The lemma is obvious from the definitions of Z and E(S). ■ 

Theorem 4.15. Take a signed graph S. Letting T range over all subsets of E, or merely 
over all for which G?(£)|T has no coloops (that is, all complements of flats of the dual bias 
matroid G _L (E) ), 

M-k) = J2\^ TC )\^\r(k+l) (4.3) 

T 

and 

^(fc) = ^MO,T c V°| T (fc + l), (4.4) 

T 

where /i is the Mobius function o/G J -(S) and 6 is the set of coloops of G(E). 

Proof. The polytope and the arrangement are transverse because intersects P°. 

Since <£> S | T = if G(E)\T has a coloop by Lemma 14.131 the two ranges of summation are 
equivalent. For S G by Lemma f4. 141 we know that 

p n s = [-1, if n ([NuiH(s|E(5) c )] x {o} B(5) ). 

Take T = E(S) C . Then 

pns=([-i,ifnz') x {o} TC 

where Z' is the real cycle space of E|T. Its Ehrhart polynomial equals (p^ T (k + 1). 

Now the result follows from Lemma 14.121 Theorem 12.21 Equation (J3.2)) . and Theorem 

EHd). ■ 
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A reminder: to apply the theorem to a graph T, take E = +I\ 

Thus, the strict integral flow polynomial can be expressed in terms of the weak polynomial 
and invariants of G _L (S). If the weak polynomial were as simple as in the case of colorings, 
where it is a monomial Section 5], we would have a nice formula for the number of 
nowhere-zero /c-flows. But such is not the case. 

It may be helpful to list some characterizations of the edge sets that support nowhere- zero 
integral flows. 

Proposition 4.16. For T C E := E(E), the following properties are equivalent. 

(i) G(E)|T has no coloops. 

(ii) E|T has a totally cyclic orientation. 
(Hi) E|T has a nowhere-zero integral flow. 

(iv) E|T has a nowhere-zero real flow. 

(v) T c is closed in the dual bias matroid G ?_L (E). 

(vi) T = E(S) C for some flat S e L(CK E ) Z ). 

Proof, (i) (ii) for graphs is Robbins' theorem jTTj. We could prove it for signed graphs, 

but the simplest approach is via oriented matroids. We know that the number of totally 
cyclic reorientations of an orientation of a matroid M is £m(0, 2) f^lE] and that this equals 
if and only if M has a coloop. Apply that to the natural orientation of G(E). 

(i) (hi) by Lemma [4.131 

(hi) ==>- (iv) is trivial. 

(iv) (i) by the proofs of j3J Lemmas 2.4 and 2.5], which amount to saying that any 
flow on E with values in an abelian group where 2a = a = must be zero on every 
coloop. Here the group is the additive group of IR. 

(v) (i): By matroid duality the complements of the closed sets in G _L (E) are the 
edge sets that do not contain a coloop of G(E). 

(v) -<=>• (vi): This is Lemma [4.121 ■ 
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